clear all 

use "data\in\elbileiere2019_maalepunkt2019_wattmete2018.dta", clear 
keep Navn Merke* Modell* Måle* Modår* mergedon 
keep if inlist(mergedon, "Navn", "name1Adresse1")
drop mergedon 
g n = _n 
reshape long Merke Modell Målepunktnr Modår, i(n)
drop n _j 
drop if Målepunktnr == . 
destring Modår, g(year)
merge m:1 Merke Modell Modår using "data\in\elbiltyper_Modell_Modår_Antall_Batterikapasitet_kWh.dta", keepusing(Batterikapasitet) nogen 
replace Batterikapasitet = 16 if Batterikapasitet == . // assume missing is small battery 
su Batterikapasitet, detail 
g byte abovemedianbattery = (Batterikapasitet > `r(p50)')
g byte EVlarge = abovemedianbattery 
drop year 
tab Merke EVlarge, m 
keep Målepunktnr EVlarge Batterikapasitet
rename Målepunktnr maalepktnr
gsort maalepktnr -Batterikapasitet // keep larges EV battery in household 
keep if maalepktnr != maalepktnr[_n-1]
tempfile EVlarge
save `EVlarge', replace 

use "data\out\pre_patterns_sample.dta", clear 

fmerge m:1 maalepktnr using "data\out\estsample_maalepktnr.dta", keep(match) nogen 

keep if elcar == 1 
merge m:1 maalepktnr using `EVlarge', keep(master match) nogen 
tab  EVlarge elcar, m 
replace EVlarge = 0 if EVlarge == .
replace EVlarge = 0 if elcar == 0  


gunique maalepktnr if elcar == 0
gunique maalepktnr if elcar == 1 & EVlarge == 0 
gunique maalepktnr if elcar == 1 & EVlarge == 1 
gunique maalepktnr
keep if (year == 2019 & month < 12)
g byte IHD = inlist(stratgroups,2,5,8,11) == 1

g consumption_all = forbruk
g consumption_0_15 = forbruk if inrange(hour,0,15) 
g consumption_16_21 = forbruk if inrange(hour,16,21) 
g consumption_22_23 = forbruk if inrange(hour,22,23)

fmerge m:1 maalepktnr using "data\in\maalepktnr_grunnkretsincome.dta", keep(master match) nogenerate
g income = Snittprhusstand / 9.5377 


forval i = 0/1 {
    estpost su consumption* IHD income if EVlarge == `i'
	est store EV`i'
}
esttab EV1 EV0  using "output\table_B1a.txt", cells("mean (fmt(2)) sd (fmt(2) par)") varwidth(17) mtitles(EVlarge EVsmall) tex replace

gunique maalepktnr if elcar == 1 & EVlarge == 0 
gunique maalepktnr if elcar == 1 & EVlarge == 1 



*** Small and large EV 
use "data\out\pre_patterns_sample.dta", clear 

fmerge m:1 maalepktnr using "data\out\estsample_maalepktnr.dta", keep(match) nogen 

gunique maalepktnr if elcar == 0
gunique maalepktnr if elcar == 1
gunique maalepktnr
keep if (year == 2019 & month < 12)
g byte IHD = inlist(stratgroups,2,5,8,11) == 1

g consumption_all = forbruk
g consumption_0_15 = forbruk if inrange(hour,0,15) 
g consumption_16_21 = forbruk if inrange(hour,16,21) 
g consumption_22_23 = forbruk if inrange(hour,22,23)

fmerge m:1 maalepktnr using "data\in\maalepktnr_grunnkretsincome.dta", keep(master match) nogenerate
g income = Snittprhusstand / 9.5377 

forval i = 0/1 {
    estpost su consumption* IHD income if elcar == `i'
	est store EV`i'
}
esttab EV0 EV1 using "output\table_B1b.txt", cells("mean (fmt(2)) sd (fmt(2) par)") varwidth(17) mtitles tex replace

gunique maalepktnr if elcar == 0
gunique maalepktnr if elcar == 1


